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ABSTRACT 



^ ■ A convenient way to calculate iV-particle quantum partition functions is by con- 

fining the particles in a weak harmonic potential instead of using a finite box or peri- 
odic boundary conditions. There is, however, a slightly different connection between 
partition functions and thermodynamic quantities with such volume regularization. 
This is made explicit, and its origin explained to be due to the system having a 
space-varying density in an external potential. Beyond perturbation theory there is 
a potential pitfall with the method, which is pointed out. 
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In order to simplify the computation of quantum partition functions one would prefer 
to utilize a geometry which allows for separation of center-of-mass motion and preserves as 
many symmetries of the model as possible. Using a finite box to confine particles does not 
meet the first requirement, and periodic boundary conditions are in conflict with the usual 
rotiational invariance of particle interactions. A rather convenient solution is to confine 
the particles in a harmonic potential, and mimick the infinite volume limit by letting this 
potential vanish. By comparing the partition function for a single particle in a finite box 
of (D-dimensional) volume V, 

Z[ v) = V A T D = V {m/2-Kl3h 2 ) D/2 , with (3 = l/k B T, 
and in the harmonic oscillator potential $(r) = |mc<j 2 r 2 , 
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one finds a volume V u corresponding to a given potential strenght mu 2 : 

V u = (2tt//W) D/2 . (1) 

can be used as an effective volume in the calculation of e.g. virial coefficients, as was 
first done by Comtet et. al.[l] for the second virial coefficient of the anyon[2,3] gas. Their 
result is in agreement with an earlier calculation by Arovas et. al. [4] (who used confinement 
in a circular disk), but the analysis is much simplified. The same method has been applied 
by Blum et. al.[5] to a similar anyon model, and is also used in many other analyses of 
anyon systems[6, 7]. But there is nothing in the method which restricts its applicability to 
such systems only, — one expects it to be a convenient method for volume regularization in 
general. 

The connection between the thermodynamic quantities and partition functions is how- 
ever different in these two cases. This fact, and its origin, does not seem to have been 
explicitly commented on in the literature. To review briefly, all intrinsic thermodynamic 
quantities can in principle be derived from the fugacity expansion of the pressure 

oo 

(3p = Y,hz\ (2) 

£=1 

where z = exp (/?//), with \x the chemical potential. For instance, the density is given as 

oo 

p = [3zdp/dz = Y / tb i z e , (3) 
t=\ 

from which the virial expansion can be obtained (by inverting (3) into an expansion for 
z = z(p), and inserting the latter into (2)). I.e. 

oo 
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with 
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Likewise the free energy per volume unit / is found from the relation 

Pf = p log (z( P )) -pp = p log (&r x p)+J2 jzrr ae pe - 



With confinement in a finite box of volume V the connection between n-particle partition 
functions and the 6/s follows from the relation 



ex P ((3pV) = Zf = exp ]T / t^(Z{ 
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I.e., by comparing (4) and (2), 



Vb e = U e (z{ v \...,Z ( />). 



(4) 
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Explicitly, the first few U^s are 



Ui = Z u U 2 



IZ2 — z\ 



3Z3 — 3Z\Z 2 + Z\ 



- [4Z A - AZ X Z 3 - 2Z\ + 4Zj% - Z\ . 
However, with confinement in a harmonic potential equation (5) must be changed to 

v u tr D i*b t = u e (z[ u \..., zf\ (6) 

where D is the number of space dimensions. The reason for the additional factor i~ D l 2 
on the left hand side is that we are dealing with a system of spatially varying density and 
pressure. The relation (4) should thus be modified to 
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exp ( / d D rPp(r)) = f) z n = exp (f^z?U e ) , 

VJ ' n=0 \i=l / 

while the relations between local pressure p(r), density p(r), and fugacity z become 

00 . 00 p 

Mr) = £ bt (ze-fW) , P (r) =^£b e (ze~^) . (8) 
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These relations are independent of the form of the potential $(r), provided the latter is 
sufficiently well-behaved. By inserting (8) into the left hand side of (7), with $(r) = 
|mwV, and performing the space integration one obtains the announced equation (6), 
with V u as defined by equation (1). The assumption behind this derivation is that the 
potential has neglectible variation on the scale of the correlation length of the system, so 
that the local equations of state are equal to those in a homogeneous bulk system of the 
same density. 

As a simple consistency check one may compute the partions functions for ideal Bose 
and Fermi gases in the two geometries. The U^s are then found to be 
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for ferm?ons res P ec tively. Here the sum runs over all single particle energy levels, with 
the energy of level v. 
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V A T D £ D l 2 (periodic boundary conditions) 



(9) 



-D 



(harmonic potential) 



Combined with equations (5) and (6) respectively one obtains the same expression for 
hi = (±1)^ +1 A^ D £~( 1+D / 2 ) in the two cases. However, one notes a difficulty with the 
method in the case of bosons below the Bose condensation point. Above Bose condensation 
the local equations of state can be written in the parametric form 



(3p(r) 
p(r) = 



r(f ) A £] 1 [°° dr r D ' 2 - 1 log (l - z e^*W e" T ) , 

J 



r(f)A£ 



drr ' 2 - 1 



e /3$(r)+r _ z 



When $ = at its minimum the fugacity z is restricted to z < 1. This leads to a maximal 
value Pmax(^) for the local density in space dimensions D > 2, and this value is smaller than 
the bulk critical density where $(r) > 0. To increase the density further a macroscopic 
number of particles must Bose condense into the ground state. These condensed particles 
will be spatially distributed as 
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p (r) — N { - 
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with N the number of condensed particles, and r c = (h/muj) 1 ^ 2 . The "coherence length" 
r c diverges as w — > + , but much slower that the linear system size, r u = V^ D = 
(2%/ flmui 2 ) 1 / 2 : Tc/r^ = (2ir/ phto) 1 ^ 2 . In the coherence region around the bottom of the 
potential the local equations of state are different from those in a homogenous bulk medium 
of the same density. This may be attributed to long range correlations in this region, 
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< p (r) p (s) ) - ( p (r) ) ( po(s) > = N ^— J 

There is an alternative, more cumbersome but also more general, method to handle 
spatially non- homogeneous systems, based on the computation of local correlation func- 
tions. It is instructive to verify that this procedure leads to results in agreement with the 
simple arguments above. In the language of second quantization the density is given by 

p(r) = ^ f (r) \P(r), 

while the local pressure can be expressed in terms of the trace of the stress tensor 1 



T j (r) = 



^-{(di&irj) (d^(r)) + (dj&irj) (d^(r)) - i^V 2 (V( r )*(r)) } , (10) 



1 With a momentum density P l {x) = — t^h |^Jd (d^) — (<9j\£d) ^ , and use of the equations of 
motion, one finds that this expression for the stress tensor leads to a force law of the expected 
form: j^pi + d-f^ = F* , with F?'(r) = -* f (r) *(r) dj$(r). 
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i.e. 
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Here \I/(r) is the second quantized Schrodinger field (assumed to obey either Bose or Fermi 
statistics), — thus p{r) and Ty(r) are quantum operators who's thermal average we want 
to compute. This information is contained in the correlation function 



G(r,s; P) = Tr 

e -/3(ei- M ) 
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where the sum over i runs over all single particle energy levels, with energies q (we assume 
tmin > /-O and normalized eigenfunctions t/v Me = [muj / '2nh sinh(£ (3hu)] c ^ 2 , and r, s are 
D-dimensional vectors (their product means scalar product). It follows that 



P(r) = ( P(r) ) p = G(r, r, 0) = ]T (±l) m ^ A/i e"^ 2 , 
where r^ 2 = (^) tanh(|.£/3ftu;). Further 
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P(r) = 2^{^ 9[ ^ G(r ' S;/5) 



-V 2 G(r,r;/3) 



(2mu/h) D 2 rj + 2(1 - D)r 
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In the limit u — > + equations (14) and (13) turn into the corresponding expressions in 
equation (8). The parameter e which determines the first correction to this limit expression 
for the pressure is seen to be 

d 

e = £f3huj, with £[3p = z—p = p. 

oz 

In conclusion, the harmonic oscillator regularization is a convenient method to impose 
particle confinement in a finite volume. Taking into account the modifications due to the 
resulting spatial nonhomogeneity of the system all bulk thermodynamic quantities may in 
principle be calculated, at least perturbatively in a fugacity or density expansion. Non- 
perturbatively the method may break down in regions with long range correlations. 
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